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ABSTRACT 

Lines from HC3N and isotopic substituted species in ground and vibrationally 
excited states produce crowded millimeter and submillimeter wave spectra in 
the C-rich protoplanetary nebula CRL618. The complete sequence of HC3N 
rotational lines from J=9-8 to J=30-29 has been observed with the IRAM 30m 
telescope toward this object. Lines from a total of 15 different vibrational states 
(including the fundamental), with energies up to 1100 cm~^, have been detected 
for the main HC3N isotopomer. In addition, the CSO telescope has been used 
to complement this study in the range J=31-30 to J=39-38, with detections 
in five of these states, all of them below 700 cm^^. Only the rotational lines 
of HC3N, in its ground vibrational state, display evidence of the well known 
CRL618 high velocity velocity outflow. Vibrationally excited HC3N rotational 
lines exhibit P-Cygni profiles at 3 mm, evolving to pure emission lineshapes at 
shorter wavelengths. This evolution of the line profile shows little dependence 
on the vibrational state from which they rotational lines arise. The absorption 
features are formed against the continuum emission, which has been successfully 
characterized in this work due to the large frequency coverage. The fluxes range 
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from 1.75 to 3.4 Jy in the frequency range 90 to 240 GHz. These values translate 
to an effective continuum source with a size between 0.22" and 0.27" and effective 
temperature at 200 GHz ranging from 3900 to 6400 K with spectral index between 
-1.15 and -1.12. We have made an effort to simultaneously fit a representative 
set of observed HC3N lines through a model for an expanding shell around the 
central star and its associated HII region assuming that LTE prevails for HC3N. 
The simulations show that the inner slowly expanding envelope has expansion 
and turbulence velocities of ~ 5-18 kms~^ and ~ 3.5 kms~^ respectively, and that 
it is possibly elongated. Its inclination with respect to the hue of sight has also 
been explored. The HC3N column density in front of the the continuum source 
has been determined by comparing the output of an array of models to the data. 
The best fits are obtained for column densities in the range 2.0-3.5-10^^ cm~^, 
consistent with previous estimates from ISO data, and T^ in the range 250 to 
275 K, in very good agreement with estimates made from the same ISO data. 

Subject headings: stars: post-AGB-stars: carbon-rich - stars: circumstellar mat- 
ter - stars: individual: CRL618 - ISM: molecules - radio hues: stars 



1. Introduction 

CRL 618 is probably the best example of a C-rich protoplanetary nebula (PPNe) with a 
thick molecular envelope (Bujarrabal et al. 1988) surrounding a BO star and an ultracompact 
HII region from which UV radiation impinges on the envelope. Its distance of 1.7 kpc makes 
this object one of the best suited sources for studying the evolutionary stages from the 
Asymptotic Giant Branch (AGB) to planetary nebulae (Herpin ct al. 2002). The brightening 
of the HII region in the 1970s (Kwok & Feldman 1981; Martin-Pintado et al. 1993) and 
the discovery of high velocity molecular winds (HVMW) with velocities up to 200 kms~^ 
(Cernicharo et al. 1989) illustrate the rapid evolution of the central star and its influence on 
the circumstellar ejected material. From observations of molecular gas at arcsecs resolution, 
the high velocity molecular outflow and the slowly expanding envelope have been resolved 
and precisely located around the ultracompact HII region (Neri et al. 1992, Cox et al., 2003, 
Sanchez-Contreras & Sahai, 2004). 

The recent discovery of the polyynes HC4H and HCgH and of benzene {CqRq), the 
flrst aromatic molecule detected outside the solar system, in CRL618 (Cernicharo et al., 
2001a,b), emphasizes the fact that the copious mass ejection, toward the end of the AGB 
phase, and the related shock and UV driven chemistry, make the C-rich nebulae of this type 
very efficient factories of organic molecules (see detailed chemical models for this object in 
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Cernicharo 2004). Small hydrocarbons and pure carbon chains are formed and they can 
in turn be the small nuclei from which larger C-rich molecules, such as the Poli-Aromatic 
Hydrocarbons (PAHs) can be built. These molecules are proposed as responsible for the 
observed emission labeled as Unidentified Infrared Bands (UIBs). 

In order to draw the clearest picture of the chemical composition, structure and dy- 
namics of CRL618, it is necessary to perform complete line surveys in those regions of the 
electromagnetic spectrum where most molecules do emit. One of these line surveys has been 
carried out by us at millimeter wavelengths with the IRAM 30m telescope (Cernicharo et 
al., in prep.). One of the most important features seen in the survey is the presence of 
clusters of lines with a frequency separation of ~ 9 GHz, corresponding to rotational transi- 
tions in the ground and vibrationally excited states of HC3N. This molecule belongs to the 
cyanopolyynes, HC2„+iN, an important group of interstellar molecules (detected in space up 
to 2n-|-l=ll, Bell et al. 1997). Because of the relatively low-energy of their bending modes, 
they can be found in space in many vibrationally excited states. As an example Wyrowski 
et al. (1999) detected HC3N in eleven vibrationally excited states with energies up to 1100 
cm~^ toward the ultracompact HII region GlO.47+0.03. The small frequency spacings of 
the HC3N rotational transitions make this molecule an interesting tool to probe the physical 
conditions of molecular gas in C-rich circumstellar shells such as in CRL618. The HC3N line 
profiles in this object show the interesting behavior of evolving, with frequency, from P-Cygni 
to normal emission. This is presented and analyzed in this paper using a model, under LTE, 
for an expanding envelope of azimuthally symmetric geometry, around a continuum source. 

In section 2 we present the IRAM-30m and CSO observation procedures. The spec- 
troscopy of HC3N is briefly summarized in section 3. The results of the observations in- 
cluding discussions about the continuum emission and the velocity components of the HC3N 
lines are found in section 4. In section 5 we present the model that has been built to explore 
the physical conditions. Simulations are presented and discussed. Finally, the conclusions 
are presented in section 6. 

2. Observations 

The CRL 618 observations presented in this paper are part of a line survey, carried out 
with the IRAM-30m telescope, and complemented with the Caltech Submillimeter Observa- 
tory (Cernicharo et al., in prep). 

The IRAM-30m telescope has presently the largest effective aperture operating at the 3, 
2 and 1.3 mm atmospheric windows. The observations started in 1994 and were completed 
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in the winter of 2001-02. In order to avoid confusion between the signal and the image 
sidebands, the receivers were optimized for image sideband rejection larger than 12 dB, 
which resulted in practically no sideband contamination, except for just a few cases where 
a very intense line (CO, HCN,...) was present in the image sideband. Receivers covered the 
following frequency ranges (in GHz): 82-116, 130-184, 201-258, 240-279. The pointing and 
focus were always checked on the target source, since CRL618 exhibits a continuum of about 
0.3 K in the scale of the IRAM-30m telescope, across the surveyed frequency range (see 
figure 1). The pointing was kept within 2" accuracy. We carried out the observations using 
the wobbler-switching mode with offsets of 60" and frequencies of 1 Hz, in order to obtain 
very flat basehnes. This is very important because some spectra are so crowded with hues 
that their reduction could be difficult otherwise (not enough baseline). The backcnds used 
were two 512 MHz filter banks connected to the receivers operating below 200 GHz, and a 
512 MHz autocorrelator with channel width of 1.25 MHz connected to the others. System 
temperatures were typically 100-400 K at 3 mm, 200-600 K at 2 mm and 300-800 K at 1.3 
mm. 

The observations performed in several observing runs in 2000, 2001 and 2002 with the 
10.4 meter dish of the Caltech Submillimeter Observatory (CSO) at the summit of Mauna 
Kea (Hawaii) made use of a helium-cooled SIS receiver operating in double-sideband mode 
in the frequency range 280-360 GHz. The continuum emission from CRL618 at the CSO 
has a level of 0.05-0.06 K in the T^ scale, too low to allow pointing on the source. The 
pointing was therefore checked on the available planets and was kept within 3-4" accuracy 
(the HPBW ranged from 20" to 25"). The system temperatures varied from 500 K to 1200 
K. In order to obtain flat baselines wc used the chopping of the secondary mirror technique 
at a frequency of about 1 Hz and offset of 90". In order to discriminate whether the lines 
are from the signal or from the image sideband, the observations were made twice with a 
frequency shift of a few tens of MHz (usually 50). 



3. HC3N spectroscopy 

The rotational spectrum of cyanoacetylene, HC3N, has been investigated in vibrational 
states up to about 1750 cm^^ by Mbosei ct al. (2000). The resulting data base for this 
molecule has been provided to us by A. Fayt. This work has been recently published (Fayt 
et al. 2004). the frequencies have been checked against our millimeter wave line survey with 
a remarkable agreement in the detected vibrational states. 

There are three doubly degenerate bending modes of HC3N, called 2/5, z/g, and z/y in 
order of decreasing energy (223, 499 and 663 cm~^). When considering the rotation, their 
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degeneracy is broken due to the direction of the rotation axis. This is known as ^-doubhng. 
1^4 is the lowest energy stretching mode (880.6 cm~^). The other 3 vibrational (stretching) 
modes have energies above 2000 cm~^ and are not seen in CRL618. In this work, we have 
been able to detect HC3N pure rotational lines in the ground and 14 different vibrationally 
excited states. Detections have been also obtained in a number of isotopomers, but they will 
be discussed in a separate work. 

The lines are labeled by the rotational quantum numbers {Jup^Jiow), fhe vibrational 
quantum numbers {v4^,v^,Vq,Vj) and the ^-doubling parameters: [ £5 £7 ]■ Table 1 provides 
the observed line parameters (velocity, width and integrated area) of a pair of rotational 
transitions (J=12-ll and J=26-25) for all detected vibrationally excited states. The J=12- 
11 lines always display an absorption feature centered at ~ -27.5 kms~^ accompanied by 
weaker emission (undetected in some cases) centered around the LSR velocity of the source. 
In the J=26-25 lines the absorption is barely detected and the emission is always much 
stronger. The selected rotational numbers are representative of the two different types of 
line profiles observed in HC3N and in other molecules toward CRL618. 



4. Results 

Of the complete frequency survey, as described in section 2, the data set used for the 
analysis presented in this paper consists of: 

• The continuum measured in each setting across the frequency ranges covered: 80-116.5, 
132-177 and 204-276 GHz. 

• Those spectra containing HC3N rotational lines in a subset of vibrational states, se- 
lected to sample the range of the vibrational energies among those detected. The 
selected vibrational states are: (0000), (0001), (0002), (0010), (0100), (0003), (0011), 
(1000), (0012), and (0102) which vibrational energies are: 0, 223, 446, 499, 663, 669, 
722, 881, 945, and 1109 cm-^ respectively. 



4.1. Continuum emission 

The central star has an effective temperature of 30000 K and it is sourrounded by an 
HII region and its PDR. Continuum emission from the HII region and the dust affect the 
population of the molecular levels and has to be taken into account in modelling the source. 
The properties of the measured continuum will be discussed in this section. 
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The level of continuum emission is important to explain the features seen in the line 
profiles of many molecules in CRL618. Although the wobbler-switching method used at the 
30m telescope for this work should provide in principle a precise continuum level for each 
spectrum directly, the continuum levels of different spectra never match perfectly. Moreover, 
earlier work has shown evidence of variability over several years (see below). Therefore, we 
have decided to use all the pointing data to retrieve the continuum by fitting gaussians to 
the pointing scans and storing the peak of the fitted gaussian (in scale) and the central 
frequency. The obvious bad scans and those where the pointing errors were larger than 
2" have been ignored. The results are shown in figure l(open circles). Tmb and fiuxes 
have been derived from a linear fit to those data (actually a horizontal line), scaled by the 
corresponding beam and forward efficiencies of the 30m telescope. In principle, this analysis 
could be affected by the presence of lines. However, the total flux from the lines represents 
less than 3-5 % of to the continuum in most frequency settings. Wc have just ignored the 
few frequency settings with very strong lines such as the lowest rotational transitions of CO 
and HCN. In total flux our measurements give 1.75 Jy at 90 GHz, 2.4 Jy at 150 GHz and 
3.4 Jy at 240 GHz. 

Several continuum measurements of CRL618 at different frequencies exist. In an early 
interferometric mapping carried out at cm wavelengths with the VLA it was shown that the 
compact radio source had a size of 0.4" xO.l" (Kwok & Bignell 1984). Martin-Pintado et al. 
(1993, 1995) using the VLA at 23 GHz found an elliptical source with size 0.40" x0.12" and 
integrated flux density of 189±20 mJy (peak continuum fiux: 16.7 mJy beam~^) in 1993 
and 250±20 mJy in 1995. Thorwirth et al. (2003) found 26.3±3 mJy at 4488.48 MHZ and 
750±75 mJy at 40 GHz with a size of 0.34" x0.16" (peak continuum flux: 140 mJy bcam~^). 
At 1.3 mm, Walmsley ct al. (1991) found fluxes varying from 1.0 to 3.3 Jy between 1987 
and 1990. If this variability is real, it could explain part of the scattering in the data points 
presented in flgurc 1 (obtained over several winters). Earlier work (Kwok & Feldman 1981, 
Spergel et al. 1983) also present the source as variable at centimeter wavelengths. Other 
measured fiuxes found in the literature include Shibata et al. (1993), 1.48 Jy at 115 GHz, 
and Martin-Pintado et al. (1998) who measured the following at 22.3, 87.6, 147.0 and 231.9 
GHz: 0.24±0.02, l.liO.l, 1.8±0.2, 2.0±0.3 Jy, values about 25-32 % smaller than our own 
results (see figure 1) that represent an average over at least 4 years. 

Prom the line survey currently under way in the 280-360 GHz range, the measured 
continuum emission averages 0.06-0.07 K in the T^ scale (~ 0.1 K in Tmb)- This result 
agrees with the trend observed at the higher end of the IRAM 30m survey (275 GHz) given 
the smaller size of the CSO antenna (10.4 m). 
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4.2. HC3N ground vibrational state 

The profiles of HC3N rotational lines within the ground vibrational state (see figure 2) 
show evidence of the high velocity molecular wind (HVMW), also seen in lines of CO and 
other molecules. In HC3N, the HVMW component of the line profile extends some 200 kms~^ 
in some cases (see figure 3). In addition to this one, there is another component associated to 
a slowly expanding molecular envelope around the ultracompact HII region. This component 
has a maximum half power width of about 18 kms~^ and is approximately centered at the 
LSR velocity of CRL618 (-21.3 kms~^). Two different absorption dips are present (see figure 
2). The first one appears at -37 kms~^ in all the lines at least up to J=37-36. It turns 
out that this feature is related to the [ 1 0]E component of the rotational transitions of 
the (0100) state. The other absorption feature at -55 kms~^ is due to expanding gas in the 
line of sight of the continuum source and is also found in low-J HCN (Neri et al. 1992) 
and CO (Cernicharo et al. 1989, Herpin et al. 2002). On the red side of the fine, between 
approximately -8 and 7 kms~^, there is evidence of emission that Neri et al. 1992 attributed 
to an unresolved structure located to the North-East of the HII region. 

4.3. Vibrationally excited HC3N 

4.3.1. The (0001) vibrational state 

The HC3N rotational fines in the (0001) vibrational level, at least below J^j, ~ 20, 
still trace gas at higher velocities than that seen in rotational lines corresponding to higher 
energy vibrational levels (see below). For example, the J=15-14 transition in (0001) shows 
an absorption dip that clearly starts below -50 kms~^, more than 30 kms~^ apart from the 
LSR velocity of the source (see figure 3). However, for J„p above 20 this colder and faster 
outfiow is not significant in the level population anymore. 

4.3.2. Other vibrational states 

The HC3N lines in the rest of the detected vibrationally excited states display patterns 
with many common characteristics. The lines appear dominantly in emission at 1.3 mm 
(at the IRAM-30m radiotelescope) and at shorter wavelengths (CSO observations at 0.850 
mm) with a half power width ranging from 5 to 16 kms"'^, but mostly around 8-10 kms~^. 
However, the majority of lines in the 2 and 3 mm windows show P-Cygni profiles, with the 
absorption part being systematically deeper at lower frequencies. This absorption is centered 
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in the range -26 to -31 kms~^, as expected for an expanding envelope with velocities around 
5-10 kms~^ surrounding the central continuum source at the LSR velocity of CRL618. The 
absorption is less deep and shifts to more negative velocities as the emission part gets stronger 
with increasing J. At frequencies well below 80 GHz, almost all lines, from other species also, 
appear strictly in absorption. The velocity of the absorption feature in HC3N rotational lines, 
with Jup < 7, in different vibrationally excited states, observed with the Effelsberg 100m 
radiotelescope (Wyrowski et al. 2003) is around -26 or -27 kms~^. Thorwirth et al. 2003 
also report HCN AJ=0, ^-doubling transitions purely in absorption at the same velocity, 
as do Martin-Pintado et al. (1995) for NH3 (3,3). The lowest transitions we have observed 
(J=9-8, J=10-9) in some highly excited vibrational states such as (0100) and (1000) do not 
show emission and the absorption feature is also centered around -27.5 kms~^. 



5. Discussion 

In this work we assume that the bulk of vibrationally excited HC3N emission comes 
from a compact envelope expanding at about 10 kms~^ around the central continuum source. 
The availability of rotational ladders in so many different vibrationally excited HC3N states 
should provide a quite precise picture of the size ratio of the slowly expanding gas region 
with respect to the continuum source, the shape of the envelope, the temperature profile, 
velocity field and HC3N density distribution. All this is explored in sections 5.1 and 5.2. 



5.1. Model 

In order to reproduce the observed rotational HC3N ladders we have built a model, for 
LTE conditions, of an expanding envelope around a continuum source. Because of blending 
with the selected HC3N lines for the analysis, we have also introduced HC5N with an abun- 
dance ratio of 1/3 respect to IIC3N (Cernicharo et al. 1987; Guelin and Gernicharo 1991), 
and isotopic species of IIG3N using a ^^G/^^G ratio of 15. 

The lines are formed as follows: First, the central continuum source is responsible for 
the observed continuum level. The gas situated between the continuum source and the Earth 
absorbs that continuum and also emits. The gas in regions outside the line of sight to the 
continuum source contributes by emission only. The continuum source is considered opaque 
to the emission of the gas behind it. The spectrum of the continuum source and its size 
ratio with respect to the gas region, turn out to be key parameters to fit the simulations to 
the observed HG3N rotational fines. The absorption takes place primarily at the terminal 
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velocity of the expanding envelope in front of the continuum source. The emission can help 
to "fill" the absorption dip depending on the characteristics of the velocity field and the 
turbulence velocity. 

The simplest model considers a spherical expanding envelope of gas surrounding a cen- 
tral continuum source (also spherical). The envelope is divided into different shells where 
the total HC3N number density and the temperature vary according to a power law, and 
radial expansion and turbulent velocities are linearly interpolated between input values at 
the boundaries. Such a simple model can only approximately reproduce the main observa- 
tional characteristics of the HC3N rotational ladders in CRL618. A more complex model has 
been developed to account for an elongated envelope with non-radial velocity components. 
We consider an external angular size of the (spherical) slowly expanding envelope 6g and 
then we truncate it so that, in the truncation direction, the size is: 9d—0g x r^, being the 
truncation parameter (see figure 4) . Other geometric details of this model are the following: 
A first coordinate system (XYZ) is defined with the Z axis running in the direction where 
the sphere is truncated, and X, Y are orthogonal to each other and to Z. The velocity field 
will be a function of [(x^-|-y^) 2 ,z] (azimuthal symmetry in the XY plane), and the temper- 
ature and gas densities will only be a function of r=(x^-|-y^-|-z^)5 (radial symmetry). The 
second coordinate system has the axis X' in common with X just described, and the remain- 
ing two (Y',Z') are rotated respect to (Y,Z) by an angle i so that Z' points to the Earth 
(hue of sight). The "edge-on" case therefore corresponds to i=90°. The continuum source 
is spherical with an angular radius 9c. The value of 9c has to be in the range 0.15" -0.4" in 
order to be compatible with data presented in section 4.1. We use as input parameter its 
ratio Tc respect to the external angular size of the slowly expanding envelope, rc=6c/9g. The 
effective temperature of the continuum source would depend on that size in order to match 
the observed continuum. In general, the equivalent size 6c would be a function of frequency. 
However, a constant size and a frequency dependent Tc, characterized by a spectral index s 
(see following equation), fit the data very nicely (see figure 1 and table 2). 

Tcont{i^) = Tcont{2QQ GHz)[^//200]^ (1) 

In the elongated envelope case, the velocity field cannot be exactly radial: v — Vrf + 
Vxyfxy (vxy 7^ if r^ <l). The extra velocity component in the XY plane makes the velocity 
pattern look as in figure 4. This kind of velocity field has the advantage of having gas 
with the same velocity projection in regions in and out the column between the continuum 
source and the Earth. Some extra gas is thus allowed to contribute filling the absorption dip 
(not only the turbulence velocity helps to do this, see figure 5). The velocity field remains 
unchanged by rotation respect to the Y or Z axis. Therefore, in the i=0'',90'' cases ("edge- 
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on" or "orthogonal" envelope cases) , only the radiative transfer results as a function of the 
impact parameter p (defined as an angle) have to be used in the numerical convolution 
with the antenna beam pattern. If i takes any intermediate value, the numerical integration 
becomes dramatically more time consuming (10-100 times) because the azimuthal symmetry 
respect to the XY' plane is broken and the convolution has to be performed in both p and 
GxY' (azimuthal angle around the Z' or line-of-sight axis, see figure 4). 

As we have seen, the model has an important number of parameters. Some of them 
can be fixed with some confidence, based on data such as the continuum fiux (see figure 
1) and the estimated distance to the source. The other parameters, velocity field, kinetic 
temperature, shape of the slowly expanding envelope, ratio between its size and that of 
the continuum source, and HC3N column density, can be constrained by running a number 
of models. It is quite straightforward to find the best values for all parameters, with the 
exception of the kinetic temperature and HC3N column density. This is discussed below. 

5.2. Constraining the physical parameters 

The first step in the simulation process has concentrated on the continuum source. We 
have considered different sizes for it (constant with frequency) and we have calculated its 
effective temperature at a reference frequency (200 GHz) and the spectral index, s, to match 
the observed continuum (figure 1). Based on the references given in section 4.1, the range 
of sizes we have explored is 0.15" to 0.40". The results are given in table 2. For each set of 
values, the fit to our continuum data is very good, typically as shown on figure 1. As the 
external size of the gas region we are modeling is around 1.5" and the ratio of this size to 
the continuum source size is a key parameter for reproducing the observed line profiles, the 
preferred values among the five given alternatives are highlighted in bold face. 

As the number of physical parameters and lines observed arc both large, a general 
numerical fit to constrain all physical parameters simultaneously is not possible. Rather, it 
is necessary to run a large number of models in several steps in order to find a reasonable 
match to all the observed HC3N rotational ladders. As this is very time consuming, the 
simulations have been restricted to the (0002), (0010), (0100), (0003), (0011), (1000), (0012), 
and (0102) vibrational states, as explained in section 4. 

In the search for the best solution we also need to define a parameter to evaluate the 
quality of the fit. Ideally, this parameter has to give equal weight to the different lines, 
independently of the frequency resolution. We thus defined an overall x ^ follows: 
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Y,. \Yd-Ym\^ 



(2) 



where Ni is the number of suitable hues in figure 6, Yd{j) are the data (total fiux divided 
by the contiuum fiux), Ym{j) are the model results, Nch is the number of channels where 
Ym{j) 7^ 1.0, i.e. we do not consider channels where the signal is just the continuum level, 
or lines from species other that HC3N or HC5N. In addition, if a line from other chemical 
species overlaps with those from HC3N or HC5N, the corresponding spectrum is discarded. 
Note also that in figure 6 not all the £-type lines are shown in each case due to the sometimes 
large velocity separation. Only the data shown in the figure are used to calculate %. Blending 
between different vibrational states or within the same vibrational state has been treated in 
our simulations. 

The two parameters that wc have tried to determine by looking at % in a grid of models 
are the temperature of the gas and the HC3N column density. In order to illustrate how 
the best solution (6) is found we show in figure 7 the values of parameter x as a function 
of these two parameters. The overall set of parameters characterizing the slowly expanding 
molecular envelope is then shown on table 3. 

The kinetic temperature found here (250-275 K) is in very good agreement with the 
value derived from ISO IR data by Cernicharo et al. 2001a,b. The average column density 
of absorbing HC3N in front of the continuum source is found to be in the range 2.0-3.5-10^^ 
cm~^. HC3N column densities in front of the continuum source where determined to be 
5-10^^ cm-2 by Cernicharo et al. (2001a) from ISO observations of HC3N (0100) and (0001) 
bending modes in absorption around 14 /im. The difference between both estimates could be 
related to two facts: First, the ISO data lacks the spectral resolution needed to fully resolve 
the individual lines used for the estimate. Second, Cernicharo et al. 2001a assumed the same 
filling factor for the continuum source and the absorbing gas. However, dust emission at the 
wavelength of the (0100) mode of HC3N, and in particular at that of the (0001) mode could 
contribute from a region larger than the zone where HC3N is produced. The IR ISO data 
could give the same HC3N column density as that derived here from the millimeter lines, if 
the continuum level to be considered for the former calculation (originated from behind the 
gas) is only about 30-50 % of the total observed continuum fiux at 14 //m. 



6. 



Conclusions 



In this paper we have shown that toward the protoplanetary nebula CRL 618: 
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• HC3N is detected in its ground and other vibrational states with energies up to 1100 
cm~^. The J up range surveyed is from 9 to 30 (IRAM 30m telescope, 15 vibrational 
states detected) and 31 to 39 (CSO telescope, detections in 5, possibly 6 vibrational 
states). The observed frequencies match very well those provided by the model of Fayt 
et al. (2004). 

• The line profiles of the HC3N ground vibrational state show evidence of the HVMW 
already found in the strong lines of CO, HCN and HCO+ with the same velocity 
span (~200 kms^^). Other vibrational states do not show this component in their 
rotational lines, except for the (0001) state that shows it marginally in the lowest 
J lines. Absorption features in the rotational lines of the ground vibrational state 
appear at similar velocities to those of features already seen in CO and HCN. The only 
exception is the one centered at ~-37 kms~^ that is due to blending with (0100) [1 
0]E lines. 

• The pure rotational transitions in vibrationally excited states of HC3N show an evolu- 
tion in their line profiles with increasing J going from almost pure absorption centered 
at around -27 or -28 kms"^, to P-Cygni profiles and, finally pure emission centered 
quite close to the LSR velocity of the source. The intermediate-J P-Cygni profiles 
have the emission part also centered quite close to the LSR velocity, while the mini- 
mum of the absorption part shifts to more negative velocities and becomes less deep 
as J increases. The exact value of J„p at which the absorption disappears changes only 
slightly with the vibrational state: It happens at slightly higher values as the energy of 
the vibrational state increases. For example, no absorption is observed in the J=24-23 
transition of the (0001) vibrational state, whereas some absorption is still seen in the 
J=26-25 line of the (1000) one. Some rotational hues within the (0001) state still trace 
larger velocities than those in the other vibrationally excited states. 

• An model for the slowly expanding envelope, under LTE, built to explore the physical 
parameters that explain the extensive set of HC3N rotational lines detected (~300, 
although the fit has been done to a subset of 116) has shown that: 

1. The size of the inner continuum source is ~0.27" with an effective temperature 
of ~3900 K at 200 GHz and a spectral index of -1.12. 

2. The external size of this slowly expanding envelope is ~ 9—1.5". 

3. The expansion velocity field has a radial component ranging from 5 to 12 kms"-*^ 
with a possible extra azimuthal component reaching 6 kms~^ at 6'= 1.5". 

4. The turbulence velocity is ~3.5 kms~^. 

5. The temperature of the envelope is in the range 250-275 K. 
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6. The HC3N column density in front of the continuum source is in the range 2.0- 
3.5-10^'^ cm-2 (the best fit is obtained for 2.8-10^^ cm'^). 
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Table 1: Observed line parameters of the J=12-ll and J=26-25 transitions in all vibrationally 
excited states of HC3N (sec text for notation) detected toward CRL618 (vi_2: position in 
kms~^, Avi 2 width in kms~^, Ai^2 area in K-km^^; 1 refers to the emission part and 2 to 
the absorption part of the line profile). In addition to the vibrational states show here, the 
(0111) at more than 1200 cm~^ is tentatively detected. 



V. state type of 
[E cm~^] t transition 


J=12-ll 

V (GHz) vi Avi Ai V2 Av2 A2 


J=26-25 
vi V (GHz) Avi Ai 


(0001) [ 1]F 
[223] [ 1]E 


109.598825 -20.9 14.9 3.6 -31.6 22.3 -3.3 

A r\f\ A A c\ "i 11 1 1 -,i /i-~ii-~i-i/-i\r/~i-i/-ii 1 

109.442011 blended with (0010) [ 1 0]F 


237.432241 -20.8 18.7 7.2 
237.093375 


(0002) [ 0]E 
[426] [ 2]E 
[ 2]F 


109.865960 blended with [ 2]E,F 
109.870288 blended with [002]F,[000]E 
109.865960 blended with [002]E,[000]F 


237.968847 -23.7 16.0 5.9 
238.053887 -21.6 11.3 3.5 
238.010133 -21.5 10.6 3.2 


(0010) [ 1 0]E 
[499] [ 1 0]F 


109.352748 -21.3 7.4 0.47 -29.8 7.3 -0.90 
109.438718 blended with (0010) [ 1 0]F 


236.900342 
237.086494 


(0100) [ 1 0]E 
[663] [ 1 0]F 


109.183008 blended with ground state of HC3N 
109.244206 -21.3 11.0 0.35 -28.4 6.6 -0.66 


236.529566 blended with grot 
236.661402 -20.5 8.4 0.4^ 


(0003) [ 1]E 
[669] [003] E+F 
[ 1]F 


110.050809 -21.3 9.0 0.37 -28.7 6.9 -0.59 
110.211793 -21.3 4.1 0.02 -29.5 7.7 -0.79 
110.366468 -21.3 14.6 0.7 -28.3 6.2 -0.62 


238.401147 -21.2 7.47 1.0 
238.771331 -21.3 12.0 1.8 
239.082313 -20.9 11.5 1.4 


(0011) [ 1-1]E 
[721] [ 1-1]F 
[ 1 1]F 
[ 1 1]E 


109.736687 ~3 MHz away from [ 1-1]F 
109.740137 ~3 MHz away from [ 1-1]E 
109.749784 ~3 MHz away from [ 1 1]E 
109.752621 ~3 MHz away from [ 1 1]F 


237.681318 -22.1 9.2 1.5 
237.713612 blended with 
237.783921 -21.8 9.9 1.9 
237.814870 -22.9 8.6 1.4 


(1000) [ 0]E 

[880] 


109.023292 -21.3 12.7 0.22 -27.9 6.2 -0.40 


236.184052 -23.6 8.8 0.3, 


(0101) [ 1 0-l]F 
[886] [ 1 0-l]E 
[ 1 1]F 
[ 1 1]E 


109.558099 -21.3 6.0 0.07 -28.0 4.5 -0.25 
109.549318 -21.3 7.5 0.16 -27.5 3.9 -0.18 
109.564687 -21.3 7.5 0.16 -27.5 5.5 -0.29 
109.553091 -27.5 3.9 -0.18 


237.336037 -21.3 12.0 0.7, 
237.409821 blended with 
237.716883 blended with 
237.898318 -19.8 11.4 1.0^ 
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Table 1: Continued. 



V. state 
[E cm- ' 



type of 

(' Iraiisitioii 



(GHz) 



J=12-ll 



Av. 



A, 



Vl 



J=26-25 

;/ (GHz) Avi A] 



(0004) 
[892] 



0]E 

2]F 

[004]E+F 
2]E 



110.543693 

110.548972 
110.554367 
110.562578 



-21.3 10.9 0.2 -27.4 5.2 -0.35 



-21.3 7.1 8.9 
-19.9 10.5 0.6 



-28.6 4.8 -0.23 



239.370230 

239.125296 
239.510881 
238.962431 



-21.1 



(0012) 
[944] 



1 0]E 
0-1 2]F 
2]E4-F 
0-1 2]E 

1 01 F 



109.989997 
110.035645 
110.097528 
110.148765 

110.189752 



-21.3 
-21.3 



5.7 
5.0 



-27.4 
-28.2 
0.09 -27.8 
0.06 -27.7 



5.2 
4.5 
5.6 
4.7 

3.0 



-0.30 
-0.18 
-0.34 
-0.16 

-0.18 



238.254341 
238.388217 
238526.950 
238.631641 



21.7 


8.1 


0.5 


20.2 


9.3 


0.6 


22.5 


8.5 


0.5 


21.5 


8.8 


0.7 


21.2 


9.5 


0.7 



(0020) 
[998] 



0]E 
2 0]E 

2 0]F 



109.522442 

109.616120 
109.616295 



-21.3 
-21.3 



7.5 
7.2 



0.23 
0.24 



-27.4 
-29.4 

-28.8 



4.3 
7.7 

8.3 



-0.17 
-0.38 

-0.41 



237.270155 

237.468947 
23.7470747 



-21.3 

-19.1 



8.2 

7.5 



(1001) 
[1103] 



1]E 
1]F 



109.306687 
109.469407 



-27.0 
-27.0 



6.0 
6.0 



-0.17 
-0.12 



236.795297 
237.145656 



-20.0 4.5 0.1 
-21.3 9.5 0.2 



(0102) 
[1109] 



1 2]E 
1 0]F 
1 2]E 
1 2]F 
1 0]E 
-1 2]F 



109.789853 
109.847627 
109.905021 
109.868821 
109.953634 
109.989048 



blended with Ci«0 



237.803492 
237.956672 
238.164654 
237.995859 
238.178301 
238.318310 



-19.9 5.3 0.2 
blended with 

-20.3 5.6 0.2 
blended with 

-21.8 8.6 0.3 



(0005) 
[1115] 



1]E 
1]F 
5]E+F 
3]E-hF 



110.655045 
111.131162 
110.895474 
110.902009 



-21.3 
-21.3 



11.7 
15.6 



0.1 
0.26 



-21.3 2.7 0.02 



-28.2 4.8 

-27.5 2.8 

-28.2 4.8 

-28.2 5.0 



-0.13 
-0.07 
-0.13 
-0.13 



230.471871 
231.456018 
231.009985 
231.031080 



-21.3 7.9 0.2 
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Table 2: Set of models to describe the central continuum source of CRL618 according to the 
parameters defined in equation 1. The preferred is in bold face (see figure 1). 

size (") Tcont(200 GHz) spectral index s 



0.18 8600 -1.15 

0.22 6400 -1.15 

0.27 3900 -1.12 

0.32 3200 -1.12 

0.36 2300 -1.12 



Table 3: Model parameters of expanding envelope that provide the best fit to the data shown 
in figure 6. 



j)araiiieter 




units 




1.5 


arcsec 


I'd 


0.7 




i 


90 


degrees 


Tgas (at 9=9c) 


263 


K 


Tgas (at 0=9g) 


263 


K 


Vr (at 9=9c) 


5.0 


kms~^ 


Vr (at 9=9g) 


12.0 


kms~^ 


Vxy (at 9^9c) 


0.0 


kms~^ 


Vxy (at 9^9 g) 


6.0 


kms~^ 


vturb (at 9=9c) 


3.5 


kms~^ 


vturb (at 9=9g) 


3.5 


kms~^ 


[HC3N] (at 9^9c) 


154 


cm~^ 


d 


-1.8 




HC3N col. density 


2.0-3.510^^ 


cm~^ 



at p<Tc 
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Fig. 1. — Open circles: Millimeter- wave continuum of CRL618 observed by the IRAM 30m 
telescope (T^ scale, from pointing observations on CRL618 in continuum mode). Horizontal 
dashed line: Fit of these data to a horizontal line. Inclined dashed line: Fit translated into 
Tmb- Crosses: Simulated continuum Tjv^s by our model (see sections 5.1 and 5.2). 
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Fig. 2. — HC3N pure rotational lines toward the proto-planetary nebula CRL618 in 
ground vibrational state. The thick dotted line marks the LSR velocity of the source, 
thin dotted lines mark interesting features in the line profile discussed in the text. 
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Fig. 3. — Velocity components of HC3N lines toward CRL618. In the ground vibrational state 
(top two spectra) we see broad wings associated to the HVMW discovered by Cernicharo 
et al. 1989 and a narrower emission due to the slowly expanding gas envelope around the 
HII region. A blending with HC3N (0100) causes the absorption feature at ~-37 kms~^. 
The absorption centered at -55 kms~^ is also seen in other abundant molecules such as HCN 
and CO. In vibrationally excited states (bottom two spectra) the HVMW component is 
not visible. At low-J the absorption of the slowly expanding gas in front of the continuum 
source is seen centered at ~-27 kms~^. However, the absorption dissapears at J > 20 for the 
particular case of the (0001) state (the J number where this happens slightly depends on the 
vibrational state). Still, rotational lines within the (0001) state trace larger velocities than 
those in higher vibrational states, although it is not the HVMW. 
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Fig. 4. — Model of CRL618 including the description of the expanding elongated structure 
considered to account for the observed HC3N line profiles. See text for details about the 
different parameters of the model. 
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Fig. 5. — Fit of observed HC3N rotational line toward CRL618 and separate contribution 
from gas in and out of the line of sight to the continuum source. The physical parameters 
correspond to those listed in table 3. 
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Fig. 6. — Observed HC3N pure rotational lines toward the proto-planetary nebula CRL618 
in several vibrationally excited states from 446 to 1109 cm~^ compared to the results of our 
model. In a few cases the observation is not shown because the corresDonding HCqN hne is 
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Fig. 7. — Value of x (defined in equation 2) for the grid of models run to find the temperature 
of the CRL618 slowly expanding envelope and the HC3N column density that best fits the 
data shown on figure 6. 



